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Abstract 

The spectral efficiency of a representative uplink (with appropriate normaUzation) in interference-limited, spatially- 
distributed wireless networks with hexagonal cells and linear Minimum-Mean-Square-Error estimation is found to 
converge to an asymptotic limit as the numbers of base-station antennas N and wireless nodes go to infinity. A 
C ' simple approximation for the mean spectral efficiency is also found for systems with both hexagonal and random 

cells when transmit power budgets are large. It is found that for large N in the interference-limited regime, the mean 
spectral efficiency is primarily a function of the ratio of the product of N and the ratio of base-station to wireless 
node density, indicating that it is possible to scale such networks by linearly increasing the product of the number 
of base-station antennas and base-station density with wireless node density. This work is useful for designers of 
wireless systems with high inter-cell interference because it provides simple expressions for spectral efficiency as a 
function of tangible system parameters like base-station and wireless node densities, and number of antennas. These 
results were derived combining infinite random matrix theory and stochastic geometry. 
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^ ■ I. Introduction 



It is increasingly common for multiple wireless networks to be within interfering distance of each other in urban 
Q environments today due to proliferation of systems such as city-wide wireless internet access, pico-cells for mobile 
^ telephony, and wireless local-area networks. Antenna arrays at base stations can significantly increase data rates in 
such systems. It is thus important to study the spectral efficiencies (b/s/Hz) of wireless links with multiple antennas 
. ^ environments that have high base-station and wireless node densities. In such systems the densities of nodes (both 
^ in-and out-of-cell) and their distribution in space are important factors as they influence inter-node distances and 
^ hence signal and interference strength, which directly impact the Signal-to-Interference-Plus-Noise-Ratio (SINR), 
- - - spectral efficiency and ultimately data rates. 

Most works on wireless networks with multi-antenna base-stations do not explicitly consider out-of-cell interfer- 
ence as resulting from spatially distributed in-band interferers. For example, Dai and Poor [2] used random matrix 
techniques similar to those used here to obtain asymptotic expressions for the spectral efficiency in multi-cellular 
environments where multiple base-stations cooperate to jointly decode signals. Recently, [3] analyzed the capacity 
region of multi-user MIMO channels with correlated channels in the asymptotic regime when the number of nodes 
is constant but the number of transmitter and receiver antennas go to infinity. However, neither of these works 
models path-loss as a function of distance and thus do not capture the distribution of interference resulting from 
spatially distributed nodes. Monte-carlo simulations were used in [4] and [5] to analyze small, spatially-distributed 
multi-antenna cellular networks. 
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are not necessarily endorsed by the United States Government. This research was supported in part by the National Science Foundation 
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Cellular networks with single-antenna base-stations and spatially distributed nodes have been analyzed in works 
such as [6], [7], [8], [9], [10], [11] using stochastic geometry to model the spatial distribution of nodes. An extensive 
set of mathematical techniques to analyze such networks can be found in [12] and [13]. 

Stochastic geometry has also been used to study ad-hoc wireless networks with both multi and single antenna 
nodes in works such as [14]. Ad-hoc wireless networks with spatially distributed multi-antenna nodes have been 
analyzed in [15] which derived an asymptotic expression for the spectral efficiency of multi-antenna links in ad- 
hoc wireless networks as a function of path-loss exponent, link length and the ratio of the number of receiver 
antennas to node-density, and [16] which found that it is possible to linearly scale the network spectral efficiency 
density by linearly increasing the density of transmitting users with the the number of receiver antennas using a 
partial-zero-forcing receiver. More recently, [17] and [18] have found exact expressions for the CDF of the SINR 
of ad-hoc wireless networks in Rayleigh fading with MMSE receivers. The key difference between these works 
and this paper is that this paper explicitly models link lengths resulting from a cellular architecture and uses power 
control that is based on the distance of nodes from their respective base-stations rather than assuming constant link 
lengths and transmit powers. Please see [19] for an comprehensive survey of works utilizing stochastic geometry 
in both cellular and ad-hoc wireless networks. 

Here, we show that with appropriate normalization, the spectral efficiency of a representative uplink for a network 
with hexagonal cells and base-stations with N antennas using the linear Minimum-Mean-Square-Error (MMSE) 
receiver converges with probability 1 and derive an asymptotic expression for the mean spectral efficiency. We 
consider interference from spatially distributed in-cell and out-of-cell wireless nodes that have single antennas 
and transmit simultaneously in the same channel using distance-dependent power control. Additionally, we find 
an approximation for the mean spectral efficiency when base station locations form a planar homogenous Poisson 
Point Process (PPP) with area density pt. We assume that signal power decays with distance r as r~", with the 
path-loss-exponent a > 2. The wireless node density is p^, and the mean per-link spectral efficiency is expressed 
as a function of A^, p^, d, and a. While the exact cumulative distribution function (CDF) of the spectral efficiency 
would be ideal, computation of this quantity is notoriously difficult for cellular systems with power-control as the 
transmit powers of nodes depend on their location on the plane. The asymptotic analysis we use here helps alleviate 
this problem and the expressions for the mean spectral efficiency are useful in understanding the large-scale behavior 
of such networks such as the rate of spectral efficiency growth with the number of antennas. 

In the process of deriving the results for cellular networks, we show in Section |lll] that the spectral efficiency 
(with normalization) of a fixed-length, multi-antenna link in ad-hoc wireless networks with nodes transmitting at 
random, independent and identically distributed (IID) power levels with a continuous Probability-Density-Function 
(PDF) that may be mildly dependent on node position (as defined in ([2])), converges with probability 1 to a limiting 
value. This strengthens an earlier result we reported in [15] which restricted the transmit powers to a finite number 
of discrete values, independent of node position. We combine this result and the PDF of transmit powers that arise 
from hexagonal and Poisson cells, and a simple power control algorithm to find the mean-spectral efficiency of 
interference-limited links in such systems. 

These results were derived combining stochastic geometry and infinite random matrix theory, specifically the 
techniques presented by Bai and Silverstein in [20]. We validated the results for finite systems using Monte Carlo 
simulations that were also used to characterize the spectral efficiency for a given outage probability. 

II. System Models 

In this section, we describe the two system models for the main results of this paper that follow in Sections JII] 
and|IV]. We first present the system model for a network with wireless nodes transmitting at random power levels, 
followed by a model for a wireless network with base-stations which we shall also call tethered nodes. 

A. Purely Wireless Network 

Consider a planar wireless network with n wireless transmitters located at random IID points in a circle of radius 
R such that 



(1) 



We shall consider the spectral efficiency of a representative link between a receiver placed at the center of the 
circular network and an additional transmitting node (node-1) located at a distance ri as illustrated in Figure |2] 
Note that the base-stations at hexagonal lattice sites represented by the solid circles should be ignored for this part. 
The remaining n transmitting nodes numbered 2, 3, • • • , n + 1 are interferers to this link. 

Suppose that each transmitting node uses a random, IID power level with PDF fp{p) and all nodes transmit 
simultaneously in the same frequency band. Let Pi equal the transmit power of node-i and the received power due 
to node-i at a distance from the representative receiver at the origin pi = GtPir^'^, where q > 2 is the path-loss 
exponent which we assume is a rational number for technical reasons. We allow Pi to be dependent on as long 
as the following holds as R,n, N ^ oo such that n/N = c > and ([T]) always hold. 
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Q requires the transmit power of a node randomly distributed with uniform probability in the circular network be 
asymptotically independent of its normalized distance from the origin as the quantity in the brackets in Q equals 
the probability that ri/^/N exceeds {Pi/x)~. Both the hexagonal-cell and Poisson cell models have this property 
as shown in Appendix |b1 

We further assume that the representative receiver has an array of antenna elements and each wireless node 
has an isotropic antenna. We assume frequency-flat fading with circularly symmetric complex Gaussian channel 
coefficients between all pairs of antennas. Let y be an A-element vector of sampled received signals at the N 
antennas of the representative receiver. Let the (n+ 1) x 1 vector x contain the transmitted signals from node-1 and 
the n interferers with entries normalized to unit power, and the x 1 vector w contain zero-mean, IID complex 
Gaussian noise terms of variance a"^ denoted by CA/'(0, fi^) . This system can be represented by the following 
equation: 

y = Hx + w (3) 

where the Ax (n + 1) matrix H is the channel matrix whose ij-th entry is the channel coefficient between 
transmitting node j and antenna element i of the receiver. Let the i-th column of H equal ^/pigi where gj has 
IID CJ\f{0, 1) entries. Thus, gj captures the Rayleigh fading and pi captures the combined transmit power and path 
loss associated with node-i. We assume that the base stations use linear Minimum-Mean-Square -Error (MMSE) 
estimators, and the transmitting nodes use Gaussian codebooks resulting in Gaussian interference. Note that the 
linear MMSE receiver is the optimal linear receiver for Gaussian signals as it maximizes the Signal-to-Interference- 
plus-Noise-Ratio (SINR) (e.g. see [21]) which maximizes the spectral efficiency. 

For technical reasons, we assume a noise power that is a function of A^ as follows: 

^2 = (A^-f ) (4) 

where is a constant. This assumption enables the asymptotic analysis of the SINR as A — oo. Without this 
assumption, as A — )• oo, the MMSE receiver will suppress interference to levels comparable to the thermal noise 
which results in a system that is not interference limited. Defining the thermal noise power as in (H) makes the 
thermal noise power increase at the correct rate such that the system remains interference-limited. Note that while 
this assumption is necessary to make precise statements about the SINR, it is does not impact the approximations 
to the SINR and spectral efficiency in interference-limited regime which is the focus of this work. 

B. Base-station Architecture 

Consider a plane divided into cells with base stations at arbitrary locations with effective area density pt. Each 
cell is associated with one base station and is the region of the plane that is closer to that base station than any 
other. Thus, each cell is the Voronoi cell associated with its base station. We assume that all cells are bounded. 
Suppose the network model from the previous section is overlaid on this network of base-stations such that one of 
the base stations is the representative receiver. Figure shows one such case where base stations are at hexagonal 
lattice sites separated by distance d. Assume that the i-th wireless node transmits data to its nearest base station 
located at a distance ru away with power Pi where 
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min ( g-rf„PM ) . (5) 



Thus, the i-th wireless node tries to achieve a target received power (relative to path-loss) of pt at its nearest base 
station, subject to a maximum power constraint Pm- 

For a set of base stations locations, the link lengths and hence transmit powers are independent random variables 
as they depend solely on the locations of the wireless nodes, which are independent by assumption. Hence, results 
derived using the assumptions of the previous section can be applied to this network model with the transmit-power 
PDF fp{p) arising from the base-station locations as long as holds. 



HI. Main Results For Receivers At Arbitrary Locations 

In this section, we present results for the spectral efficiency of the representative link between node-1 and the 
representative receiver at the origin of the network described in Section III- A I To avoid degenerate expressions 
in the derivation, we define a normalized SINR of the representative link at the output of the MMSE receiver: 
(3^ = Ar-"/2siNR. The SINR can then be found by multiplying (3^ by iV"/^. Using this definition, we introduce 
the following theorem proved in the Appendix |A] 

Theorem 1: Consider the network model from Section III-AI As the number of interferers n — oo, the number 
of antennas N ^ oo, and the outer radius of the network i? — )• cx) such that c = n/N > and = are 
constants, then f3j\[ — )• /3 with probability 1 where /3 is the unique, non-negative real solution to the following 
equation: 



E[P^/'^]0 



2/a 







— CSC 




a 





oo 



fp{x)x''dx dr 



0rr^ P{- ... 

TTa = 2 

2Gtp^TTP; = 2p^7rrf 

where b = (^^) ^ , E[P'^] is the expected value of the transmit power of the wireless nodes raised to ^, and Pi 
is the power of the representative transmitter. 

Since all nodes use Gaussian codebooks the residual interference at the output of the MMSE receiver is Gaussian, 
so we use the Shannon formula to estimate the per-Unk spectral efficiency of the system: 

C = log2(l + SINR) = log2(l + 7V"/2/3) . 

Since the log function is continuous, as n, N ^ oo in the manner described in Theorem 1, the following expression 
holds with probability 1 (e.g. see [22]): 

C-log2(iV-/2)^log2(/3). (7) 

Hence, with appropriate normalization, the spectral efficiency converges to an asymptote with probability 1 as 
N ^ oo. Following the steps of the derivation in Appendix E of [23] which is essentially an application of the 
dominated convergence theorem (e.g. see [22]), the deviation of the mean spectral efficiency from is asymptotic 
value approaches zero, i.e. 



£;[C]-log2(l + iV"/'/3) 



0. (8) 



Hence, log2 (l + TV"/^^) is a good approximation for the mean spectral efficiency for large N. 

From Theorem [T] the value of the limiting normalized SINR /3 is unclear as it is given implicitly by the solution 
to To obtain a more meaningful expression for /3 and the spectral efficiency, we can simplify Equation (|6]l 
using Lemma 3 of [23] which indicates that as c — )• oo, 6 — ?• and the second term on the LHS of Q vanishes. 
Thus, when the ratio of the number of antennas at the representative receiver to the number of interferers is high 
(i.e. large c), Q can be written as 0: 



'This approximation requires the solution of (|6) to be a continuous function of /3. If a is rational, ^ can be raised to a sufficiently high 
power resulting in a polynomial equation in /? with real coefficients that are known to have continuous roots. 



Additionally, in the interference-Limited regime, we assume that cr^ is sufficiently small that the second term on 
the LHS of (|9l) is dominated by the first. Writing = [^sin (^)] ^ for convenience, neglecting the second 
term on the LHS of substituting the definition of and rearranging terms yields the following approximation 
for the SINR when is large: 



SINR « PiGc 
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(10) 



We have made several approximations in deriving (ITO^ . The validity of these approximations for reasonable 
values of A^, n, and cr^ are verified in simulations presented in Section IIV-CI Results of more extensive simulations 
can be found in [24]. 

With Gaussian codebooks, the mean spectral efficiency can be approximated as: 



E[C(ri,Pi)] «log 1 + 



N 



(11) 



Suppose that the maximum distance between any transmitting node and its desired receiver < ( 



.We 



call this the sufficient-power case since every wireless node can satisfy the target received power pt at its desired 
receiver. For the cellular model described in the next section, this corresponds to the base-station separation being 
sufficiently small that the target received power is attained by each wireless node. Substituting (|5]l into ([TTI l. 
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which is a function of the second moment of the link-lengths. 



IV. Hexagonal Cells 

Suppose that the base stations are located at hexagonal lattice sites on the plane separated by distance d which 
results in hexagonal cells as illustrated in Figure [T] The following lemma statistically characterizes the link-lengths 
for this model. 

Lemma 1: The PDF fx{x), Cumulative-Distribution-Function (CDF) Fx{x), and /c-th moment of the link length 
X between a randomly located wireless node and its closest base station in a hexagonal-cellular system with minimum 
base station separation d are the following: 
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Proof: Consider Fig. [U which illustrates a portion of a wireless network with hexagonal cells. Each wireless 
node in the network falls on some random point in an equilateral triangle formed by the three base stations closest 



to it, and foms a link with the base station at the closest vertex of that triangle as illustrated in Fig. [T] Thus, the 
link-lengths are statistically equivalent to the distance between a randomly selected point in an equilateral triangle 
to the closest vertex of that triangle. The CDF, PDF and k-th moments of the distance between a random point in 
an equilateral triangle to the closest vertex are known [25], and are precisely the formulae in Lemma [T] Note that 
the PDF of link-lengths associated with a single hexagonal cell which equals ([T3] ). has been given without proof 
before in [26]. ■ 

A. Sufficient Transmit Power 

If d < -75 ( ) ° , all wireless nodes have sufficient transmit power to meet the target received power pt at 



V3 \ Pt 

their base-stations. From ([121 1. the spectral efficiency depends on the second moment of link-lengths given by ([15] ) 
with k = 2: 

Substituting into ([T2l) yields the following expression for the mean uplink spectral efficiency of interference-limited, 
hexagonal-cell systems with a large number of base sation antennas: 

E[C\^\ogAl + G^(—^, VV (16) 



B. Insufficient Transmit Power 



If d > (^^^^ j ° ' transmit power budget is insufficient for all nodes to meet the target received power at 
their base stations which results in some wireless nodes transmitting at full power. In this case, 

^[p2/a] (which IS 

required to find the mean spectral efficiency using (ITTi ) is given by the following lemma which can be proved by 
direct computation using Lemma [T] 



Lemma 2: If Pm < ^ (f) , 



E[p-.]=P^j-^i-l] P/,. (17) 
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Lemma |2] substituted into (|lll) yields the mean spectral efficiency for a link of length ri . Averaged over the PDF 
of link-lengths arising from hexagonal cells, the mean spectral efficiency is: 

E[C] = I log, I 1 + G^^x^ ( -— 1 I Ux)dx 

y Gt \E[Pc]7rp^x^ ^ 

+ 1/ . ^ log2 1 + G^Pm \ 1 I fx{x)dx 
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where Fx{x) and fx{x) are from Lemma 1, and £'[P^/°] is from Lemma |2] We were not able to integrate the 
second term on the RHS of ( [T9l ) and thus use numerical integration to compute E[C] for this case. 

C. Monte Carlo Simulations 

To verify the asymptotic results of the previous section, we simulated network topologies with base stations 
at hexagonal lattice sites, and wireless nodes distributed randomly on a large circular network on the plane. We 
simulated each configuration 5000 times. To reduce edge effects, we evaluate uplink spectral efficiencies in the 
center-most cell using the Shannon formula. 

For each trial, we placed 4000 wireless nodes randomly in circular networks with radii selected to meet target 
wireless node densities of 10^^, 10^^, and 10~^ nodes m^^. The circular network was overlayed on a hexagonal 
grid of base stations which extends beyond the edge of the circular network of wireless nodes. The base stations 
were spaced such that their densities were 20%, 10%, 5% and 2.5% of the wireless node density. 

The thermal noise power was fixed at 10~^^VF, equivalent to an antenna temperature of ~ 300 K for 200kHz 
bandwidth. The target received power at the base stations pt, was such that the received SNR = SOdS. We 
used a high value for pt to ensure that the system is interference-limited for low densities of wireless nodes. We 
simulated systems with both unlimited transmit powers (to simulate the sufficient-power case) and powers limited 
to Pm = 200ml^. 

The channel coefficient between the antenna of wireless node i and antenna j of the representative base station 

was modeled as \J Gtr^'^gij, where a = A, Gt = 10~^ m^. Pi is the transmit power of the wireless node (a function 
of the distance between the wireless node and its nearest base station) and Qij are IID CAA(0, 1) random variables 
which represents the narrow-band Rayleigh fading channel. 

1) Sufficient Transmit Powers: Figure |3] illustrates the mean uplink spectral efficiency for wireless node densities 
of Pw = W^^ and pw = 10^^ nodes m~^, and unlimited transmit powers per node versus the number of antennas 
at the representative base station. The square and asterisk markers represent wireless node densities of 10~^, and 
10~^ nodes m^^, respectively and the solid lines represent the asymptotic mean spectral efficiency from ([T6] l. 

Note that the asterisk and square markers coincide indicating that the absolute density of wireless nodes does not 
effect the mean spectral efficiency, and it is the relative density of wireless to base stations that matters. Furthermore, 
it is clear that the asymptotic approximation ([T9] l holds when N is sufficiently large. For instance, when the base 
station density is 20% of the wireless node density, the asymptotic and simulated mean spectral efficiency differ 
by less than 10% when N > 10. For lower densities of base stations, the convergence is slower, e.g. when the 
base station density is 5% of the wireless node density, the difference between the simulated and asymptotic mean 
spectral efficiency drops below 10% only when N > 37. 

We analyzed the outage spectral efficiencies from the simulated data, where spectral efficiency with outage 
probability Pq means that a fraction 1 — Pq of the links in our simulations achieved that spectral efficiency or 
greater. Figure |4] illustrates the outage spectral efficiencies vs. number of receive antennas at the representative base 
station for p^ = 10^^ nodes with 5%, 25% and 50% outage probabilities. 

Note that the intersection of the line with the circular markers and the 1 bs^^Hz^^ mark in Figure |4] occurs 
approximately at = 14 indicating that it is possible for 95% of links to achieve 1 bs~^Hz~^ with > 14 
antennas at the base stations when the base station density is 10% of the density of transmitting wireless nodes. 
In real systems, the number of nodes transmitting at any time is far smaller than the total number of nodes in the 
network. Suppose that at any one time, 10% of nodes are actively transmitting in the network. Figure 51 indicates 
that with a base station density equaling 1% of total wireless node density (including inactive ones), it is possible 
for 95% of links to achieve 1 hs~^Hz^^ with 14 antennas at each base station. 

2) Insufficient Transmit Power: Figures |5] and |6] illustrate the mean spectral efficiency vs. number of receive 
antennas for = 10~^ and p^j = 10~^ respectively, with Pm = 200mVF. The markers represent the simulated 
mean spectral efficiencies for different relative densities of base stations to wireless nodes. The solid lines are the 
predicted asymptotic mean spectral efficiencies obtained by numerically evaluating equation ([T9l ). 

It is clear from Figures |5] and |6] that the asymptotic approximation ( fT9l ) holds when N is sufficiently large. 
In Figure |5l the simulated and asymptotic mean spectral efficiencies agree to within 5% for A'^ > 2 for all the 
base-station densities considered. In Figure |6] however, for base-station densities that are 5% of the wireless node 
density of 10^^ nodes the simulated and asymptotic spectral efficiencies differ by less than 13% only when 



there are 13 or more antenna elements at the receiver. For base station densities that are 20% of the wireless node 
density, the simulated and asymptotic spectral efficiencies agree to within 13% when > 3. 

At low wireless node densities, the simulated spectral efficiencies converge more rapidly (compared to high 
densities) to the asymptotic values because a large fraction of nodes transmit at the 200 mW power limit. The 
empirical distribution function (e.d.f.) of interference powers at the representative receiver thus converges more 
rapidly to its asymptotic value. The rate of convergence of the e.d.f. of interference powers controls the rate of 
convergence of the eigenvalues of the spatial interference covariance matrix GPG^ (see Appendix |A] and Section 
3 of [27]) which affects the convergence rates of the normalized SINR and spectral efficiency. 

Figure |7] shows the outage and mean spectral efficiencies for = 10~^ (solid lines) and pw = 10"'^ (dashed 
lines) nodes 10% relative density of base-stations to wireless nodes and Pm = 200m VF. Note that with 

10 antennas at the receiver, approximately 0.2 and 0.3 b/s/Hz are achievable for p^ = 10~^ and pw = 10"^ 
respectively. The discrepancy in the spectral efficiency is a result of the maximum transmit power. For p^ = 10"*^, 
a larger fraction of nodes transmit at Pm compared to p^ = 10^'^, resulting in higher Signal-to-Interference-Ratios 
(SIR) for Pw = 10^^. The higher total interference power for p^, = 10~^ is offset by increased signal powers due 
to shorter links since the relative base-station to wireless node density is fixed. 



V. Random Cells 

A. Mean Spectral Efficiency 

Suppose that instead of at hexagonal lattice sites, the base stations were located at random points on the plane 
according to a homogenous PPP with intensity pt nodes m~^. The random cells generated constitute a homogenous 
Poisson-Voronoi tessellation (PVT) of the plane where the Voronoi cell associated with each base station is the 
set of points on the plane that are closer to that base statio than any other. We shall henceforth omit the term 
homogenous, i.e. PPP and PVT shall refer to the homogenous PPP and homogenous PVT respectively. Figure |8] 
illustrates a portion of such a network where the base stations aie the circles and the cell boundaries are the solid 
lines. The square is a representative wireless node connected to its closest base station. 

The distances between wireless nodes and their closest base station (and their transmit powers) are correlated 
as they are related by the random locations of the base stations. This correlation prevents direct application of 
Theorem 1 which requires independent transmit powers. However, conditioned on a realization of the PPP, the 
transmit powers are independent as they are simply functions of the node locations which are independent by 
assumption. Thus, we first find the mean spectral efficiency of the representative link conditioned on a realization 
of the base station PPP lit, and then average over all 11^ to find the unconditional mean spectral efficiency. First, 
place the representative transmitter at an arbitrary point on the plane. Next, consider a realization of lit and let 
the closest base station to the representative transmitter be the representative receiver. We assume that lit does not 
result in any Voronoi cell of infinite area. Excluding such realizations of PVTs does not change the mean spectral 
efficiency since they are zero-probability events (e.g., see [28]). Additionally, we shift the coordinates of the system 
such that the representative receiver is at the origin, n IID interferers are then placed in a radius- ii disk centered 
at the origin. We shall consider the limit as R, N and n — )• oo as in the previous section. 

Conditioned on lit, and ri, the mean spectral efficiency is approximated by: 
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with probability 1 (see e.g. [28] for a discussion on ergodicity). 

r — I 

That is to say, E[Po lilt J which is averaged over all the cells in the PVT, is independent of the specific realization 
of the PVT with probability 1. Hence, 
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The PDF of the distance r between any wireless node and its closest base station frir) = 2-Kptre ^p*^^ , for 
r > (e.g. see [25]). Using this expression and the power control of ([5]), 

+ P^,27Tptre-^P^^ dr = {^] . (22) 
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Substituting ([5]) into (|2TI) and averaging with respect to ri, 



e Vpt ; log2 1 + —Go ' 



Gt \E[py-]7Tp, 

+ . . log2 ( 1 + PMTi^Go, { ^^^,7a]vrpJ ^ ) 2mne--''-?dn. (23) 



We were unable to find a closed form expression for the second term on the RHS of (1231) and thus use numerical 
integration to evaluate it. However, if the transmit power budget of each wireless node is large (or the density of 
base stations is high), (l22l ) is approximated by 

E[pf^^] « (24) 



E[G] « log2 I 1 + ptGo, \ 1 I . (25) 



and (1231) is approximated by 



Substituting (1241) into (1251 ). 



E[C]^logJl + Ga(^^] ], (26) 



i.e., the mean spectral efficiency primarily depends on the ratio of pt and p^ implying scale invariance in dense 
networks where the transmit power limit Pm is not significant. 

Note that while (l26l ) does not depend on the choice of pt the original equation used to derive (|26] | was based 
on the assumption that the system is interference limited which means (|26l ) is valid only when pt and p^ are 
sufficiently high that the system is interference limited. 

The scale invariance implied by (l26l ) indicates that with hexagonal cells, constant mean spectral efficiency can 
be maintained by fixing the relative density of base stations to wireless nodes. 



B. Monte Carlo Simulations 

We verified (1231 ) and (l26l ) by Monte Carlo simulations of the network topology. We placed base stations in a 
circular network of radius SR. The numbers of base stations were selected to achieve relative densities of base 
stations to wireless nodes of 20%, 10% and 5%. The network of base stations was then re-centered such that a 
base-station exists at the origin. 4000 wireless nodes were then placed in a circular network of radius R, centered 
on the network of base stations with R selected to achieve a wireless node density of 1 x 10^^ nodes m^^. This 
experiment was repeated 1000 times. For each trial, the spectral efficiency of a randomly selected link in the 
center-most cell (to reduce edge-effects) was collected and averaged. The transmit power of each wireless node 



was set according to ([5) with Pm = oo (to simulate the sufficient power case) or Pm = 200 mW. Gt = 10 ^m", 
thermal noise power of 10"^^ W, and Q = 4, were assumed. 

Figure |9] shows results of Monte Carlo simulations and the asymptotic expression given bv(|26]l for systems with 
unlimited transmit powers per node. Note that the simulations match the asymptotic results to within 10% when 

> 9 for a relative base station to wireless node density of 20%. For lower relative densities, the convergence 
is slower. For 10% relative density, the the simulations match the asymptotic expression to within 10% only when 

> 19 and only when > 38 for 5% relative density. The rate of convergence for random cells is slower 
than that for hexagonal cells because the range of transmit powers is much larger for random cells compared to 
hexagonal cells which results in slower convergence, as explained in Section IIV-CI 

Figure [10] shows simulations of systems with a 200 mW transmit power limit. The target received power pt was 
set such that the target SNR, pt/fj"^ = 30dB. For relative base station to wireless node densities of 20%, 10%, and 
5%, the simulated mean spectral efficiencies are within 10% of the asymptotic prediction when > 6, > 7 and 
A^ > 9 respectively. The convergence of the simulated mean spectral efficiencies to the asymptotic values is faster 
for systems with limited transmit power as the range of transmit powers in the network is smaller when there is a 
bound on the transmit power. 

These simulations indicate that the asymptotic expressions are useful for reasonable numbers of base-station 
antennas. 

C. The Cost of Random Cells 

For systems with limited transmit powers, we numerically evaluated and plotted equations for the spectral 
efficiency corresponding to random and hexagonal cells in Figure \TT\ where the solid and dashed lines represent 
hexagonal and random cells respectively. The transmit power budget was 200 mW and wireless node density was 
10"'^ with different relative density of base stations to wireless nodes as shown in the plot. Note that the difference 
in mean spectral efficiencies diminishes with the number of antennas. However, for high base station densities the 
mean spectral efficiency for random cells is significantly lower. For instance, with 10 antennas at the base stations 
and 20% relative density of base stations to wireless nodes, the mean spectral efficiency with hexagonal cells is 
twice that of random cells. 

When base-station density and/or transmit power budgets are high, the mean spectral efficiency given by ([T6l ) 
can be rewritten in terms of the effective base-station density as: 



When compared to the mean spectral efficiency with random cells given by ( |26l) . ( |27l) indicates that several-fold (but 
not orders of magnitude) gains in mean spectral efficiency can be achieved by regularly distributing base stations 
in planar networks compared to randomly distributing them, and furthermore, the difference diminishes with the 
number of base-station antennas. 



We have derived an asymptotic expression for the mean spectral efficiency of the uplink in wireless networks 
with multi-antenna base-stations (tethered-nodes) in networks with hexagonal cells. We assumed a power control 
algorithm for which wireless nodes try to achieve a target received power at the base stations to which they are 
connected. This power control algorithm which has also been used in [14] and related works ensures that uplink 
spectral efficiencies are close to the mean value with high probability when the number of antennas per link is 
large and the wireless nodes have high power budgets. 

If the spacing between base stations is small enough that all wireless nodes are able to achieve the target received 
signal power at their base stations (which we call the sufficient-power case), the mean spectral efficiency takes 
a simple form given by ([T6] l. Note that for a fixed ratio of base station to wireless node densities pt/pw, as 
increases, the system eventually moves to the sufficient-power case so this is an effective way of scaling the density 
of such networks. 

From ([T6l ). note that with 7 antenna elements per base station and pt/pw ~ 0.1, the mean spectral efficiency is 
approximately 1 bs^^Hz^^. If we assume that 10% of all wireless nodes are actively transmitting at any one time. 




(27) 



VI. Summary and Conclusions 



the ratio of base station to total wireless node density has to be just 1% to achieve a mean spectral efficiency of 1 
bs~^Hz~^, as given by (fT6l ). For systems with insufficient power, i.e., the base stations are far enough apart that 
some fraction of the wireless nodes will not achieve the target received power, the expression for the mean spectral 
efficiency is more complicated and has to be evaluated by numerical integration. 

We verified the accuracy of the derived expressions by Monte Carlo simulations. We also used the simulations to 
study the outage spectral efficiency, i.e., the spectral efficiency achievable with a given probability. We found that 
in the sufficient power case, with 14 antennas per base station (a reasonable number for base-stations) and single 
antennas at each wireless node, and with 10% of wireless nodes transmitting simultaneously at any one time, over 
1 bs'^Hz^^ is achievable by 95% of wireless nodes when the ratio of base station to wireless node densities is 
1%. 

We also found an expression for the mean spectral efficiency of links with base stations at random locations with 
area density pt. We find that the penalty of random cells compared to hexagonal cells diminishes with increasing N. 
At modest N we found that hexagonal cells can increase the mean spectral efficiency over random cells, several-fold 
as illustrated in Figure \TT\ 

The findings of this work are useful for designers of cellular wireless systems such as pico-cells and city-wide 
wi-fi access as they provide compact expressions for the spectral efficiency and hence data rates as a function of 
tangible system parameters such as user and base-station densities, number of base-station antennas and random 
versus regular distribution of base-stations. 



Appendix 



A. Derivation of Asymptotic Normalized SINR for General Link Lengths (Proof of Theorem I) 

Let the representative transmitter be node-1, at a distance ri from the representative receiver at the origin. 
The remaining transmitting nodes numbered 2, 3, • • • n + 1 are interferers. To find the normalized SINR, /S^ = 
A^~~SINR, we scale the interference and noise powers by N~ and find the SINR of this new system. Let pi = 
N~pi = N~ PiGtr^°^ where Pi is the transmit power of node-i, and its distance from the origin. Let the matrix 
P = diag{p2,P3, ■ ■ ■ ,Pn+i)- Since the noise power of the original system is a"^ from dU, the noise power of the 
new system is N'^l'^a^ = Na'^. By the well-known formula for SINR of MMSE estimators. 



N 



GPG^ + cj^I 



h 



(28) 



where h is an x 1 vector of channel coefficients between the representative transmitter and the antennas at 
the representative receiver, and G is an n x A matrix of IID CM{0, 1) Rayleigh fading coefficients between the 
antennas of the interferers and the representative receiver. 

By Theorem 7.1 of [20] which is a strengthening of Theorem 3.1 of [21], if the e.d.f of the received interference 
powers (i.e. P2,P3,-" iPn+i) converges with probability 1 to a limiting function H{x) as A^ and n — )• oo with 
n/N = c > 0, the SINR converges with probability one to the value of m{z) that satisfies ( |29l ) with z 

TdHir) 



zm{z) + 1 = m{z)c 



(29) 



/o 1 + Tm{z) ' 

To show convergence of the interference powers, recall that n interferers are distributed in a disk of radius R 
centered at the origin. Setting = 1 for notational convenience (it will be reintroduced in the final expressions), 
the CDF of the received power from wireless node i is 



Pr{pj < x} 



Pr{P,A-rr" 



< x} 



Pr 



X 



Pi) fp{Pi)dPi. 



(30) 




fp{P)p- dP. 



(31) 
(32) 



(1301) to dlD is from ©, and from substituting b = (^)"''^ 

By the Glivenko-Cantelli theorem, the e.d.f. of a set of IID random variables converges uniformly, with probability 
1, to its CDF. Hence, the e.d.f. of the piS converges with probability 1 to the RHS of ( l32l ). i.e. H{x) = Pr{pj < x}. 
Taking the derivative of the RHS of (l32l ) yields: 



dH{x) 2iTpu 



dx 



ca 



■E 



pi 



ca 



x/b 



fp{T)T'-dT. 



Substituting (1331) and applying Lemma 1 of [15], the RHS of ( l29l ) becomes 



mc 



TdHij) 
1 + rm 



27rp^, 



■E 



a 



pi 



mc 



m ° vr CSC 





2 " 

P- 


Jo ca 





1 + rriT 



-dr 



2-k\ 2Tr p^m 
a J a 



1 + rriT 



■dT 



fp{x)xo,dx . 



r/b 



Substituting dMll into ^ yields 



2ttp^ 

zm[z) + 1 = E 



P- 



a 

2TTpwm{z) 
a 



m[z) ° TT esc ( — 1 



l + m(z)r 



dr 



fp{x)x''dx. 



rjb 



Since /3iv — >• /3 = m(— cr^), substituting z = —a^ and /3 = m(z) into (1351 ) yields: 



V/3 + l=^i? 



/3 a vr esc 



l + /3r 



fp{x)x<^dx. 



r/b 



Re-arranging terms yields ©. 



(33) 



(34) 



(35) 



(36) 



B. Mild Dependence of Transmit Power on Position for Cellular Networks 

We show that ri/\/N and Pi are asymptotically independent for both hexagonal and Poisson Voronoi cells. 
Assume that node i is distributed with uniform probability in the radius R circular network and recall that r^, ru 
and Pi are its distance to the origin, distance to its closest base station and transmit power respectively. Let be 
a disk of radius x centered at the origin and Hj, be the union of the set of disks of radius v centered at the base 
stations. Conditioned on rj < yVN, node-i is unifomily distributed in B^^. Hence, 

r , n 1 r , Area(H,nS^) 
VTlrti<v\y=<y \ =VT^rti<v\ri<y^/N^ = ^-^^^ ^ (37) 



For hexagonal cells, as and R — )• oo, the RHS approaches Fx{v) from (lT4l i as the edge effects diminish. For 
Poisson cells, the set H^, forms a Boolean model and as — )• oo, the RHS of ( 1371 ) converges to 1 — e~f*'"^ with 
probability 1 (see Chapter 3 of [28]). Hence, ru and ri/\/N are asymptotically independent and since Pi is a 
function of ru. Pi is asymptotically independent of ri/\fN. Thus, as — )• oo, with probability approaching unity. 



Note that the RHS is simply the CDF of evaluated at \fN (^) °. Substituting c = n/N, (H]), Bayes' rule, and 
rearranging terms in the RHS of the last equation yields 
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• Base stations 

X Representative wireless node 




Fig. 1. Illustration of base stations at hexagonal lattice sites. 




Fig. 2. Illustration of wireless network with representative link and base-stations at hexagonal lattice sites. 




Number of Base Station Antennas 



Fig. 3. Mean spectral efficiency vs. number of receive antennas for p„ = 10 ^ and p„ = 10 nodes m ^ with unlimited transmit powers 
and hexagonal cells. 
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Fig. 4. Outage spectral efficiency vs. number of receive antennas for p^, = 10 ^ nodes m with unlimited transmit powers and base 
station density equaling 10% of wireless node density, with hexagonal cells. 
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Fig. 6. Mean spectral efficiency vs. number of receive antennas for — 10 ^ nodes m ^ with different relative density of base stations 
to wireless nodes and hexagonal cells. 




p 


= 50% 







p 


= 25% 







p 


= 5% 







p 


= 1% 



Number of Base Station Antennas 



10^ 



Fig. 7. Outage spectral efficiency vs. number of receive antennas for ~ 10"'' and 10^"^ nodes and base-station density equal to 
10% of wireless node density, hexagonal cells and 200mW transmit power budget. The solid lines represent p^ = 10"'^ and dashed lines 
represent puj = 10~*. The markers represent the different outage probabilities shown in the legend. 
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Fig. 8. Illustration of network with base-stations distributed according to a Poisson point process. 
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Fig. 9. Mean spectral efficiency of uplink communications with random cells and unlimited transmit powers. 
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Fig. 10. Mean spectral efficiency of uplink communications with random cells and 200mW transmit power limit per node. 
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Fig. 11. Mean spectral efficiency of the uplink with random cells and hexagonal cells and transmit power limited to 200 mW. Solid and 
dashed lines represent hexagonal and random cells respectively. 



